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Abstract 

Background: Metabolic profiling may provide insight into biologic mechanisms related to age-related increases in regional 
adiposity and insulin resistance. 

Objectives: The objectives of the current study were to characterize the association between mid-thigh intermuscular and 
subcutaneous adipose tissue (IMAT, SCAT, respectively) and, abdominal adiposity with the serum metabolite profile, to 
identify significant metabolites as further associated with the homeostasis model assessment of insulin resistance (HOMA- 
IR), and, to develop a HOMA-IR associated metabolite predictor set representative of regional adiposity, in 73 functionally- 
limited (short physical performance battery <10; SPPB) older adults (age range, 70-85 y). 

Methods: Fasting levels of 181 total metabolites, including amino acids, fatty acids and acylcarnitines were measured with 
use of an untargeted mass spectrometry-based metabolomic approach. Multivariable-adjusted linear regression was used in 
all analyses. 

Results: Thirty-two, seven and one metabolite(s) were found to be associated with IMAT, abdominal adiposity and, SCAT, 
respectively, including the amino acid glycine, which was positively associated with SCAT and, negatively associated with 
both IMAT and abdominal adiposity. Glycine and four metabolites found to be significantly associated with regional 
adiposity were additionally associated with HOMA-IR. Separate stepwise regression models identified glycine as a HOMA-IR 
associated marker of both IMAT (model R 2 = 0.51, p<0.0001) and abdominal adiposity (model R 2 = 0.41, p<0.0001). 

Conclusion: Our findings for a positive association between glycine with SCAT but, a negative association between glycine 
with IMAT and abdominal adiposity supports the hypothesis that SCAT metabolic processes are different from that found in 
other fat depots. In addition, because of the significant associations found between glycine with HOMA-IR, IMAT, SCAT and 
abdominal adiposity, our results suggest glycine as a serum biomarker of both insulin sensitivity and regional fat mass in 
functionally-limited older adults. 
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Introduction 

Insulin resistance is well documented to occur with advancing 
age [1] and, is associated with an elevated risk of developing 
cardiovascular disease [2,3,4] and mortality [5]. One explanation 
for the increased prevalence of insulin resistance found in older 
adults involves elevated levels of total body adiposity [6] . However, 
not only the total amount of adipose tissue seems to be important, 
but also its distribution in the body [7]. For example, lower 
extremity intermuscular adipose tissue (IMAT) is elevated in older 
adults [8] and, is associated with insulin resistance [9,10]. 
Similarly, age-related increases in abdominal adiposity may 
explain the increased incidence of insulin resistance found in 
older adults [11,12,13]. Interestingly, thigh IMAT and subcuta- 



neous adipose tissue (SCAT) have been reported to have opposing 
effects on insulin sensitivity [14]. Thigh SCAT has previously been 
shown to be positively associated with insulin sensitivity, when 
examined independently of the other adipose tissue depots of the 
thigh [14,15], and, elevated levels of thigh SCAT are indepen- 
dently associated with more favorable glucose and lipid profiles in 
overweight older adults [16]. Therefore, developing an improved 
understanding about mechanisms that underlie both regional 
adiposity and insulin resistance may lead to targeted interventions 
designed to reduce fat mass and to improve insulin sensitivity in 
older adults. 

One approach that can be used to elucidate mechanisms related 
to both regional adiposity and insulin resistance is use of mass 
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spectrometry-based (MS) metabolomics, which (when an untar- 
geted approach is used) aims to characterize and quantify all of the 
small molecules in a sample (blood, urine, cerebrospinal fluid, etc.), 
thereby providing an analytical description of complex biological 
media [17]. Mass spectrometry-based (MS) metabolomic signa- 
tures have proven their value in several diseases, including diabetes 
and cardiovascular disease [18], and, cancer [19]. For example, 
MS-based metabolomic profiling of amino acids, fatty acids, and 
acylcarnitines [20] has provided support for the hypothesis that 
lipid-induced insulin resistance is explained in part by overload of 
mitochondrial lipid oxidation, accumulation of incompletely 
oxidized fats, and depletion of TCA intermediates, thereby 
leading to a condition of mitochondrial stress that activates 
signaling pathways that interfere with insulin action [2 1] . To date, 
only one study has investigated the association between circulating 
biomarkers (vitamins, adipokines and hormones) with regional 
adiposity in older adults (average age 63 y) [22], and, studies 
investigating the association between metabolic intermediates with 
both regional adiposity and insulin resistance in an older 
population have yet to be performed. Therefore, the goals of the 
present study were to characterize the association between serum 
MS metabolites with thigh IMAT, SCAT, and abdominal 
adiposity, to investigate the association between significant 
metabolites with HOMA-IR, and, to develop an insulin 
resistance-associated metabolite predictor set representative of 
regional fat mass, in functionally-limited older adults. 

Methods 

Ethics Statement 

All participants signed an informed consent. Consent was 
confirmed by a witness to the process. Copies of signed consent 
were provided to all participants and originals were retained in 
study documents. All consent procedures and the study protocol 
were reviewed and approved by the Tufts University Health 
Sciences Campus Institutional Review Board. 

Participants 

Serum from 73 (40 women, 30 men) overweight older adults 
(average BMI and age, 27.0 kg/m 2 , 77.8 y) with a demonstrated 
reduction in mobility (a score of 10 or less on the Short Physical 
Performance Battery, SPPB) from the randomized, double-blind, 
controlled study of Chale et al. (2012) [23] was subjected to 
analysis. All participants were required to be sedentary, defined as 
the absence of structured exercise during the previous 6 months. 
28 of the 73 subjects in the current study were classified as insulin 
resistant (average HOMA-IR, 3.7), based on the HOMA-IR 
criteria defined by Stern et al. (2005) [24] . 

Measurement of Regional Adiposity 

Following an eight hour fast, thigh IMAT and SCAT were 
measured by computed tomography using a Siemens Somotom 
Scanner (Erlangen, Germany). Manual tracings were used to 
delineate thigh IMAT and SCAT in the range of -190 to -30 
Hounsfield units (HU) using SliceOmatic v4.2 software (Montreal, 
Canada), as described previously [25]. Average measurement, 
standard error of measurement (SEM) and the intraclass 
correlation coefficient (ICC) for IMAT were 5.9 (range 3.2-9.3), 
0.6 and 0.99 (95% CI 0.95, 1.0), respectively. Average measure- 
ment, SEM and ICC values for SCAT were 66.0 (range 26.4- 
138.4), 11.4 and 0.999 (95% CI 0.995, 1.0). To measure regional 
(arms, legs, trunk) and total fat mass, whole body dual-energy-x- 
ray absorptiometry scans were performed using a Hologic 
Discovery A densitometer. Scan acquisition and analysis were 



performed according to manufacturer guidelines, with three passes 
over the subject to acquire the full DXA image Scans were 
analyzed using Hologic QDR software version 12.3 in array mode 
(Hologic, Inc. Bedford, MA). Abdominal adiposity was determined 
via the trunk fat-periphery fat ratio, calculated by dividing the 
trunk fat mass by the sum of fat mass in the arms and legs [22]. 

Measurement of HOMA-IR 

Serum glucose was measured by the hexokinase method using a 
Beckman-Coulter AU400e chemistry system (Beckman-Coulter, 
Diagnostic Division, Brea, CA). Serum insulin was measured using 
a competitive binding radioimmunoassay (HI- 1 4K Human Insulin 
Specific Kit, Millipore, St. Charles, MO). The Homeostasis Model 
Assessment of Insulin Resistance (HOMA-IR) is commonly used 
in large epidemiological studies and in clinical practice to estimate 
insulin resistance from fasting plasma insulin and glucose 
concentrations (Matthews 1985). HOMA-IR was calculated with 
the following formula: {(fasting serum insulin (mU/L) x fasting 
plasma glucose (mmol/L))}/22.5 [26]. 

Metabolomic Analysis 

Ten mL of blood was collected under standardized conditions 
between 8 and 10 AM and, following an overnight fast. After 
collection, blood was allowed to clot for one hour at room 
temperature, and, centrifuged at 2, 1 35 xg for 10 minutes at 4°C. 
Serum was derived by removing the supernatant following 
centrifugation and was stored in 1 mL aliquots at — 80°C, prior 
to analysis. 

Metabolomic data acquisition was performed by Metabolon 
Inc. (Research Triangle Park, NC). Briefly, small molecule 
metabolites were extracted from serum and the reconstituted 
extracts were resolved using mass spectrometry platforms, 
including ultrahigh performance liquid chromatography/ tandem 
mass spectrometry (UHPLC/MS/MS) optimized for acidic and 
basic species, and, gas chromatography/mass spectrometry (GC/ 
MS), with details of this platform described extensively by Evans 
et al [27]. Of the 296 metabolites that were identified, 181 
metabolites were determined to be representative of fatty acids, 
amino acids, and acylcarnitines [20,28,29]. MS-based metabolo- 
mic profiling of amino acids, fatty acids, and acylcarnitines has 
previously been used to explain mechanisms related to lipid- 
induced insulin resistance [20,21]. To examine the association 
between metabolites previously shown to be associated with insulin 
resistance with regional adiposity and HOMA-IR, serum data for 
amino acids, fatty acids, and acylcarnitines and, their related 
metabolites were used in multivariable statistical analyses. 

Statistics 

Linear regression (SAS Enterprise Guide 4.3; SAS Institute Inc., 
Cary, NC) was used in multivariable analyses to determine the 
association between IMAT, SCAT and abdominal adiposity with 
serum metabolites, between significant metabolites with HOMA- 
IR, and, to identify an insulin resistance-associated metabolite 
predictor set representative of IMAT, SCAT and, abdominal 
adiposity. 

Separate multivariable linear regression models were performed 
for each mass spectrometry-based metabolite. All regional fat 
models included sex, age, total fat and metabolite. HOMA-IR 
models included sex, age, total fat and metabolites found to be 
significandy associated with regional adiposity. Metabolite param- 
eter estimates are based on a one-unit increase in the median 
scaled, natural log transformed metabolite data. False discovery 
rates [30] were computed with use of the q-value method [31] to 
adjust for multiple comparisons. Because conservatively chosen q- 
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values may exclude biologically relevant associations, we chose to 
use a liberal q-value cutoff of 0.30 (and, p^O.05), as reported by 
Meyers et al (2010) [32]. A q-value of 0.30 indicates that the result 
is likely to be valid 7 out of 10 times, which we suggest is 
reasonable in the setting of exploratory discovery. 

To determine an insulin resistance-associated metabolite 
predictor set representative of IMAT and abdominal adiposity, 
sex, age and total fat-adjusted stepwise linear regression was used 
on metabolites that were found to be significandy associated with 
either IMAT or abdominal adiposity and, HOMA-IR. Because 
only one metabolite (glycine) was found to be associated with thigh 
subcutaneous adipose tissue, stepwise regression was not per- 
formed for SCAT. When stepwise selection is used, metabolites 
enter and remain in the model only if they are statistically 
significant (p^0.05) after inclusion. Internal validation of the 
stepwise models was performed using a bootstrap re-sampling 
procedure of the dataset, using 200 individual bootstrap replicates. 
The pre-specified models were fitted to each, and an empiric 95% 
confidence interval (95% CI) was generated about the median of 
the 200 replicate parameter estimate values for each covariate. 

Results 

General physical characteristics, including age, percent female, 
body weight, body mass index (BMI), total fat mass, trunk fat/ 
peripheral fat, IMAT, SCAT and HOMA-IR are shown in 
Table 1. 

Metabolites Associated with IMAT 

Significant sex, age and total fat-adjusted associations between 
thirty-two metabolites with IMAT were identified, including 
eighteen amino acids and amino acid-related metabolites, eleven 
fatty acids and fatty acid-related metabolites, and, three carnitines 
(Table 2). Of the eighteen significant amino acids and amino acid- 
related metabolites, twelve metabolites were found to be positively 
associated with IMAT, including 2-hydroxyisobutryate, N-acetyl- 
threonine, urea, 4-acetamidobutaiioate, indolelactate, N-acedya- 
lanine, indoleacetate, symmetric dimethylarginine, y-glutamyl- 
threonine, N-mefhylproline, 2-hydroxy-3-methylvalerate and 
creatine. Six amino acids and amino acid-related metabolites 
were found to be negatively associated with IMAT, including 
glycine, glutamine, S-methylcysteine, serine, methionine and 
asparagine. 

All eleven fatty acids and fatty acid-related metabolites were 
found to be positively associated with IMAT, including saturated 



Table 1. Subject demographics. 




Age (y) 


77.7±3.9 


% Female 


59 


Body weight (kg) 


74.0±10.6 


BMI (kg/m 2 ) 


27.0±3.2 


Total Fat Mass (kg) 


25.6±7.2 


Trunk Fat/Peripheral Fat 


1.1 ±0.3 


IMAT (cm 2 ) 


4.7 ±2.6 


SCAT (cm 2 ) 


82.5 ±42.4 


HOMA-IR 


3.7±2.0 



Values shown represent means±SD. IMAT, intermuscular adipose tissue; SCAT, 
subcutaneous adipose tissue; homeostasis model assessment of insulin 
resistance, HOMA-IR. 
doi:1 0.1 371 /journal.pone.0084034.t001 



Table 2. Metabolites significantly associated with IMAT. 







P + SE 


p-value 


q-value 


Amino Acid Metabolism 


2-hydroxyisobutyrate 


2.5±0.6 


1E-04 


0.007 


Glycine 


-2.9±0.9 


0.001 


0.03 


Glutamine 


-5.1 ±1.7 


0.003 


0.04 


N-acetylthreonine 


2.5±0.8 


0.003 


0.04 


Urea 


1.6±0.6 


0.007 


0.07 


S-methylcysteine 


-1.7±0.6 


0.01 


0.07 


4-acetamidobutanoate 


2.2±0.9 


0.01 


0.07 


Serine 


-2.9±1.2 


0.01 


0.07 


Methionine 


-4.5±1.9 


0.02 


0.08 


Indolelactate 


1.5±0.7 


0.03 


0.10 


N-acetylalanine 


3.1 ±1.4 


0.03 


0.10 


Indoleacetate 


1.5±0.7 


0.03 


0.10 


Symmetric dimethylarginine 


2.5 ±1.2 


0.03 


0.10 


y-glutamylthreonine 


1.4±0.7 


0.04 


0.10 


N-methyl proline 


0.5±0.2 


0.05 


0.10 


2-hydroxy-3-methylvalerate 


1.2±0.6 


0.05 


0.10 


Asparagine 


-2.5±1.2 


0.05 


0.10 


Creatine 


1.1 ±0.6 


0.05 


0.10 


Fatty Acid Metabolism 


17-methylstearate 


2.1 ±0.8 


0.009 


0.07 


Dihomo-linolenate (20:3n3 or n6) 


2.4±0.9 


0.01 


0.07 


Palmitate (16:0) 


2.9±1.1 


0.01 


0.07 


Pentadecanoate (15:0) 


2.0±0.8 


0.02 


0.08 


Docosadienoate (22:2n6) 


2.0±0.8 


0.02 


0.08 


Margarate (17:0) 


1.8±0.8 


0.02 


0.08 


10-heptadecenoate (17:1n7) 


1.7±0.8 


0.03 


0.10 


Myristate (14:0) 


2.1 ±0.9 


0.03 


0.10 


Docosapentaenoate (n6 DPA; 22:5n6) 


1.4±0.6 


0.04 


0.10 


Stearidonate (18:4n3) 


1.3±0.6 


0.04 


0.10 


Caproate (6:0) 


2.0±1.0 


0.05 


0.10 


Carnitine 


Butyrylcarnitine 


1.9±0.5 


0.0006 


0.02 


3-dehydrocarnitine 


2.4±0.8 


0.003 


0.04 


Methylglutaroylcarnitine 


0.6±0.3 


0.04 


0.10 



All models were adjusted for sex, age and total fat. Metabolites are listed with 
their respective parameter estimates and standard errors (P±SE) in order of 
significance (p-value) and, with q-values. 
doi:1 0.1 371 /journal.pone.0084034.t002 



(caproate, C6:0; myristate, C14:0; pentadecanoate, C15:0; palmi- 
tate, CI 6:0, 17-methyl stearate, margarate, CI 7:0), unsaturated 
(10-heptadecenoate, 17:ln7) and polyunsaturated (stearidonate, 
18:4n3; dihomo-linolenate, 20:3n3 or n6; docosadienoate, 22:2n6; 
docosapentaenoate, n6 DPA; 22:5n6) fatty acids. 

Three carnitines were found to be significantly positively 
associated with IMAT, including butyrylcarnitine, 3-dehydrocar- 
nitine and methylglutaroylcarnitine. 

Metabolites Associated with Abdominal Adiposity 

Significant sex, age and total fat-adjusted associations between 
seven metabolites with abdominal adiposity were identified, 
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including five amino acids and amino acid-related metabolites 
and, two carnitines (Table 3). Of the five significant amino acids 
and amino acid-related metabolites, positive associations between 
2-hydroxy-3-methylvalerate and 3-hydroxyisobutyrate, and, neg- 
ative associations between glycine, glutamine and 3-methoxtyr- 
osine with abdominal adiposity were found. Among the carnitines, 
significant positive associations between isovalerylcarnitine and 3- 
dehydrocarnitine with abdominal adiposity were identified. Fatty 
acid metabolites were not found to be significantly associated with 
abdominal adiposity. Although ten additional metabolites (17- 
methylstearate, 2-hydroxybutyrate, a-hydroxyisovalerate, creatine, 
dihomo-linolenate (20:3n3 or n6), palmitate, palmitate methyl- 
ester, serylleucine, trans-4-hydroxyproline, valine) were found to 
be significant at p^0.05, q-values greater than 0.30 excluded them 
from statistical significance [data not .shown). 

Metabolites Associated with SCAT 

In contrast with the significant negative associations found 
between glycine with both IMAT and abdominal adiposity, 
glycine was found to be significantly positively associated with 
SCAT (P±SE, 35.0±9.8, p = 0.0007, q = 0.13) after adjusting for 
sex, age and total fat. Although nine additional metabolites (2- 
hydroxy-3-methylvalerate, propionylcarnitine, butyrylcarnitine, 3- 
hydoxyisobutyrate, alanine, isovalerylcarnitine, leucylalanine, 
tiglylcarnitine, trans-4-hydroxyproline) were found to be signifi- 
cant at p£0.05, q-values greater than 0.30 excluded them from 
statistical significance (data not .shown). Interestingly, multiple 
metabolites found to be positively associated with either IMAT 
or abdominal adiposity (2-hydroxy-3-methylvalerate, butyrylcar- 
nitne, isovalerylcarnine, 3-hydoxyisobutyrate) were found to be 
negatively associated (p^O.05) with SCAT. 

Metabolites Associated with Regional Adiposity and 
HOMA-IR 

Metabolites found to be associated with IMAT, abdominal 
adiposity, SCAT and, HOMA-IR are shown in Table 4. Glycine 
was found to be significantly negatively associated with HOMA-IR 
and, IMAT and abdominal adiposity. In contrast, a significant 
positive association between glycine with SCAT was identified. 2- 
hydroxy-3-methylvalerate was found to be significantly positively 
associated with HOMA-IR, IMAT and, abdominal adiposity. In 
contrast, a negative association (p£0.05) between 2-hydroxy-3- 



Table 3. Metabolites significantly associated with abdominal 
adiposity. 









P+SE 


p-value 


q-value 


Amino Acid Metabolism 


Glycine 


-0.3±0.1 


0.002 


0.15 


2-hydroxy-3-methylvalerate 


0.2 ±0.1 


0.007 


0.19 


Glutamine 


-0.5 ±0.2 


0.008 


0.19 


3-hydroxyisobutyrate 


0.3 ±0.1 


0.007 


0.19 


3-methoxytyrosine 


-0.3±0.1 


0.01 


0.24 


Carnitine 


Isovalerylcarnitine 


0.2 ±0.1 


0.001 


0.15 


3-dehydrocarnitine 


0.3 ±0.1 


0.006 


0.19 



All models were adjusted for sex, age and total fat. Metabolites are listed with 
their respective parameter estimates and standard errors (J3±SE) in order of 
significance (p-value) and, with q-values. 
doi:1 0.1 371 /journal.pone.0084034.t003 



methylvalerate with SCAT was identified, but a q-value greater 
than 0.3 (0.67) excluded it from statistical significance. Glutamine 
was found to be significantly negatively associated with HOMA- 
IR, IMAT and, abdominal adiposity. 2-hydroxyisobutryate and 
butyrylcarnitine were found to be significantly positively associated 
with HOMA-IR and IMAT. It is important to note that although 
butyrylcarnitine was found to be negatively associated with SCAT 
at p£0.05, a q-value greater than 0.3 (0.67) excluded it from 
statistical significance. 

Stepwise Regression Identifies Glycine as a HOMA-IR 
Associated Predictor of IMAT and Abdominal Adiposity 

To determine a HOMA-IR associated metabolite predictor set 
representative of IMAT and abdominal adiposity, multivariate- 
adjusted stepwise regression was used (Table 5). Although sex, age 
and total fat were found to explain 36.0% of the variability 
inherent in IMAT, the combination of glycine and 2-hydro- 
xyisobutyrate were found to explain an additional 15%, for a total 
adjusted R" of 51%. Similarly, although the combination of sex, 
age and total fat explained 33.4% of the variation inherent in 
abdominal adiposity, glycine was found to explain an additional 
8%, for a total R 2 of 41.4%. It is important to note that use of 
stepwise regression identified glycine as a negative predictor of 
both IMAT and abdominal adiposity. Furthermore, because 
glycine was the only metabolite significantly positively associated 
with thigh subcutaneous adipose tissue, stepwise regression for 
SCAT was not determined. To assess internal validity of the results 
obtained from stepwise regression bootstrapping was performed. 
200 bootstrap samples were created, and the variables that were 
found to be significant in the stepwise models were forced in. The 
median R 2 across the bootstrap replicates for the IMAT and 
abdominal adiposity models was found to be 0.55 (95% CI 0.38, 
0.68) and 0.45 (95% CI 0.27, 0.57), respectively, which is similar to 
the model performance found in the stepwise models. Bootstrap- 
based 95% CI for the metabolites identified using stepwise 
regression are shown in Table 5. 

Discussion 

The goals of the current study were to characterize the 
association between serum MS metabolites with thigh IMAT, 
SCAT, and abdominal adiposity, to investigate the association 
between significant metabolites with HOMA-IR, and, to develop a 
HOMA-IR-associated metabolite predictor set representative of 

Table 4. Metabolites associated with HOMA-IR, IMAT, 



abdominal adiposity, and, SCAT. 













Abdominal 






HOMA-IR 


IMAT 


Adiposity 


SCAT 


Glycine 


-2.4±0.8* 


-2.9±0.9* 


-0.3+0.1* 


35.0±9.8* 


2-hydroxy-3- 


1.8±0.5* 


1.2 ±0.6* 


0.2±0.1* 


-16.4±6.7* 


methylvalerate 










Glutamine 


-4.8±1.5* 


-5.1 ±1.7* 


-0.5±0.2* 


23.7±19.9 


2-hydroxyisobutyrate 


1.3±0.6* 


2.5 ±0.6* 


0.0+0.1 


-1.4±7.6 


Butyrylcarnitine 


1.5±0.5* 


1.9 ±0.5* 


0.1 ±0.1 


-15.2±6.4* 



All models were adjusted for sex, age and total fat. Metabolites are listed with 
their respective parameter estimates and standard errors (J3±SE). 
*p<0.05, q <0.030; 
*p<0.05, q>0.30. 

doi:1 0.1 371 /journal.pone.0084034.t004 
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Table 5. Stepwise linear regression on metabolites 
containing significant associations with regional adipose 
tissue and HOMA-IR. 






P±SE 


p-value 


Bootstrap p (95% CI) 


IMAT 


Intercept 


-17.3 ±4.9 


0.0008 


-16.7 (-26.3, -6.8) 


Sex 


1.2±0.5 


0.03 


1.6 (0.6, 2.8) 


Age 


0.2 ±0.1 


0.002 


0.2 (0.1, 0.3) 


Total fat 


0.2 ±0.0 


2E-08 


0.3 (0.2, 0.4) 


2-hydroxyisobutyrate 


2.0±0.9 


0.002 


2.0 (0.9, 3.0) 


Glycine 


-2.0 ±0.6 


0.03 


-1.9 (-3.5, -0.4) 


Abdominal Adiposity 


Intercept 


1.0±0.6 


0.11 


1.0 (-0.0, 2.2) 


Sex 


0.4 ±0.1 


4E-07 


0.4 (0.2, 0.5) 


Age 


-0.0 ±0.0 


0.67 


-0.0 (-0.0, 0.0) 


Total fat 


0.0 ±0.0 


0.08 


0.0 (-0.0, 0.0) 


Glycine 


-0.3±0.1 


0.002 


-0.3 (-0.5, -0.2) 



All models were adjusted for sex, age and total fat. Covariates are listed with 
parameter estimates and standard errors (P±SE), and with p-values. Bootstrap 
parameter estimates (fj) are shown with their respective 95% confidence 
intervals (CI). 

doi:1 0.1 371 /journal.pone.0084034.t005 



regional adiposity, in functionally-limited older adults. Multivar- 
iable-adjusted linear regression on serum metabolomic data 
identified thirty-two, seven and one metabolite(s) to be significantly 
associated with IMAT, abdominal adiposity and SCAT, respec- 
tively, including the amino acid glycine, which was positively 
associated with SCAT and negatively associated with both IMAT 
and abdominal adiposity. Glycine and, four metabolites found to 
be significantly associated with regional adiposity were additionally 
associated with HOMA-IR. Multivariable-adjusted stepwise 
regression identified glycine as a HOMA-IR associated predictor 
of both IMAT and abdominal adiposity. In addition, glycine was 
the only metabolite to reach statistical significance in its association 
with SCAT, thereby identifying glycine as a HOMA-IR associated 
marker of multiple adipose-containing compartments. Because 
glycine has previously been shown to reduce plasma insulin and fat 
mass in rodents [33,34], future studies aimed at testing the 
hypothesis that glycine supplementation may reduce fat mass and 
improve insulin sensitivity in glycine-deficient functionally-limited 
older adults are of interest. 

Upon first inspection of the metabolites found to be associated 
with IMAT, two main trends were evident, including positive 
associations for all eleven fatty acids and, negative associations for 
five glucogenic amino acids. Elevated levels of plasma free fatty 
acids have been shown to reduce insulin-stimulated glucose uptake 
in both young healthy subjects [35], and, in those with type II 
diabetes [36]. Furthermore, we found significant positive associ- 
ations for palmitate and dihomo-linolenate, fatty acids which are 
associated with peripheral insulin resistance in elderly men [37]. 
IMAT has been suggested to contribute to insulin resistance by 
enhancing rates of lipolysis within skeletal muscle [38], a finding 
that may explain the multiple associations found between serum 
fatty acids with IMAT. In addition, five glucogenic amino acids 
(asparagine, serine, glycine, glutamine, methionine) were found to 
be negatively associated with IMAT. Because the utilization of 
carbohydrate is somewhat impaired when insulin resistance is 
present, the oxidation of amino acids becomes an alternative 



energy source by entering the citric acid cycle at different points 
[39]. For example, asparagine is converted into oxaloacetate, 
serine and glycine are converted into pyruvate, glutamine is 
deaminated to yield ot-ketoglutarate, and methionine yields 
succinyl-CoA [40]. Although prior evidence indicating these 
amino acids as markers of IMAT has yet to be reported, serum 
levels of glutamine and methionine were reduced in subjects with 
type II diabetes, relative to subjects with normal glucose tolerance 
[41]. Similarly, methionine, glutamine, asparagine, serine and 
glycine and were each reduced in plasma from 1 1 -week old 
diabetic rats, relative to lean controls [42] . Collectively, our results 
suggest that the reduced levels of these amino acids previously 
found in association with insulin resistance and type II diabetes 
[41,42] may be related to elevated levels of IMAT. 

Separately, five amino acids and amino acid metabolites and, 
two carnitines were found to be significantly associated with 
abdominal adiposity. Interestingly, glycine and glutamine were 
negatively associated with both IMAT and abdominal adiposity, 
whereas three branched chain amino acid (BCAA)-related 
metabolites, 2-hydroxy-3-mefhylvalerate (formed during isoleucine 
degradation [43]), 3-hydroxyisobutyrate (formed during valine 
degradation [44]) and isovalerylcarnitine (the cognate acylcarni- 
tine of isovaleryl-CoA, formed during leucine degradation [45]) 
were found to be positively associated with abdominal adiposity. 
Isoleucine deamination yields 2-oxo-3-methylvalerate [43]; reduc- 
tion of this metabolite then produces 2-hydroxy-3-methylvalerate. 
Branched chain 2-keto dehydrogenase, the protein responsible for 
degradation of 2-oxo-3-methylvalerate has been shown to be 
reduced in adipose tissue from obese rodents [46]. The potential 
inability to degrade 2-oxo-3-methylvalerate may then lead to the 
accumulation of this metabolite, followed by mass action-driven 
reduction to produce 2-hydroxy-3-methylvalerate. Similarly, 
although 3-hydroxyisobutyrate and isovalerylcarnitine have yet 
to be reported as markers of abdominal adiposity, these 
metabolites are elevated in type II diabetics [47] and in obese, 
insulin resistant subjects (when compared with lean controls) [20]. 
As a possible explanation for these findings, multiple genes 
involved in BCAA catabolism have been shown to be down- 
regulated in adipose tissue from obese, when compared with non- 
obese subjects [48]. For example, isovaleryl CoA dehydrogenase 
gene expression is decreased in subcutaneous abdominal adipose 
tissue from insulin resistant, when compared with insulin sensitive 
subjects [49]. Collectively, these data suggest that mechanisms 
related to BCAA degradation in abdominal adipose tissue may be 
impaired in functionally-limited older adults. 

In contrast with the negative associations found between glycine 
with both IMAT and abdominal adiposity, glycine was found to be 
positively associated with SCAT. In addition, the negative 
associations found between 2-hydroxy-3-methylvalerate and 
butyrylcarnitine with SCAT (found just outside of statistical 
significance, p^0.05 but q ^0.030), contrasts the positive 
associations found between these metabolites with either IMAT 
or abdominal adiposity. Furthermore, thigh SCAT was negatively 
associated ((3±SE, -0.0±0.0, adjusted R 2 = 0.21, p = 0.004), 
whereas both IMAT and abdominal adiposity were positively 
associated with HOMA-IR (IMAT (3±SE, 0.3±0.1, adj. 
R 2 = 0.19, p = 0.01; abdominal adiposity (5±SE, 2.0±0.9, adj. 
R 2 = 0.17, p = 0.03), after adjusting for age, sex and total fat. In 
agreement with these data, opposite associations for SCAT and 
IMAT with insulin sensitivity in obese HIV+, and, in healthy 
postmenopausal women, respectively, have been reported [14,15]. 
As a potential explanation for the differences in insulin sensitivity 
found when comparing fat depots, adipocyte progenitor cells 
(APC) isolated from SCAT have intrinsic differences in gene 
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expression, differentiation properties, and responses to environ- 
mental/genetic factors when compared with APC from visceral 
adipose tissue [50]. In addition, transplantation of SCAT from 
donor mice into visceral fat regions of recipient mice decreased 
body weight, total fat mass, glucose and insulin levels, evidence 
that suggests that SCAT produces substances that can act 
systemically to improve glucose metabolism [51]. Collectively, 
these data suggests that mechanisms related to glycine metabolism 
may underlie the metabolic differences inherent in thigh SCAT, 
when compared with thigh IMAT or abdominal adiposity. 

Metabolites found to be significantly associated with regional 
adiposity were then further examined for their association with 
HOMA-IR, and, negative associations for both glycine and 
glutamine with HOMA-IR were identified. Glycine and glutamine 
have been reported to be negatively associated with HOMA-IR in 
subjects from the Framingham Heart Study [52]. Although 
concentrations of glycine in plasma and serum are reduced in 
insulin resistant, obese human subjects [20,42,53], our results are 
the first to more specifically associate glycine and glutamine with 
both HOMA-IR and, with IMAT, abdominal adiposity or SCAT. 
Furthermore, positive associations for 2-hydroxy-3-methylvalerate, 
2-hydroxyisobutyrate and butyrylcarnitine with HOMA-IR were 
identified. 2-hydroxy-3-methylvalerate is a product of isoleucine 
degradation [43]. Isoleucine has been previously reported to be 
among an amino acid combination that is predictive of future 
diabetes risk [54]. Although isoleucine was not found to be 
associated with HOMA-IR or regional adiposity in the current 
study, our results additionally suggest the isoleucine degradation 
product 2-hydroxy-3-methylvalerate as an insulin resistance- 
associated marker of thigh IMAT and abdominal adiposity. 
Separately, increased urinary excretion of 2-hydroxyisobutyrate, 
derived from the microbial degradation of dietary proteins [55] 
has been reported in type 2 diabetic mice [56] and in obese 
humans [57]. Similarly, plasma levels of butyrylcarnitine were 
elevated in a rodent diabetes model, relative to controls [58], and 
have been shown to be elevated when mitochondrial fatty acid 
oxidation is impaired [59]. Collectively, these results suggest that 
mechanisms related to isoleucine degradation, gut microbial 
metabolism, and mitochondrial dysfunction may underlie both 
insulin resistance and regional adiposity, in functionally-limited 
older adults. 

To determine an insulin resistance-associated metabolite 
predictor set representative of regional adiposity, we used 
multivariable-adjusted stepwise regression and found glycine 
(and 2-hydroxyisobutyrate, with IMAT) to be associated with 
both IMAT and abdominal adiposity, in separate models. 
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